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Abstract 

A toy model for magnetic extraction of energy from black hole (BH) accretion 
disk is discussed by considering the restriction of the screw instability to the mag- 
netic field configuration. Three mechanisms of extracting energy magnetically are 
involved. (1) The Blandford-Znajek (BZ) process is related to the open magnetic 
field lines connecting the BH with the astrophysical load; (2) the magnetic coupling 
(MC) process is related to the closed magnetic field lines connecting the BH with its 
surrounding disk; and (3) a new scenario (henceforth the DL process) for extract- 
ing rotational energy from the disk is related to the open field lines connecting the 
disk with the astrophysical load. The expressions for the electromagnetic powers 
and torques are derived by using the equivalent circuits corresponding to the above 
energy mechanisms. It turns out that the DL power is comparable with the BZ and 
MC powers as the BH spin approaches unity. The radiation from a quasi-steady 
thin disk is discussed in detail by applying the conservation laws of mass, energy 
and angular momentum to the regions corresponding to the MC and DL processes. 
In addition, the poloidal currents and the current densities in BH magnetosphere 
are calculated by using the equivalent circuits. 
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1 INTRODUCTION 

As is well known, the Blandford-Znajek (BZ) mechanism has been re- 
garded as a reasonable process for powering the radio jets in AGNs (Blandford 
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& Znajek 1977; Rees 1984). Recently the BZ mechanism has been used as a 
central engine for powering gamma-ray bursts (GRBs), where rotating energy 
of a stellar black hole (BH) with magnetic field of 10^^ gauss is extracted along 
the magnetic field lines supported by a magnetized accretion disk (Lee et al. 
2000; Wang et al. 2002a). According to the BZ theory the power extracted 
from the BH arises from the rotation of the open magnetic field lines relative 
to the BH horizon, resulting in the electromotive force (EMF) in an equivalent 
circuit (MacDonald and Thorne 1982, hereafter MT82). The BZ power can be 
effectively regarded as the power of the current dissipated on the astrophysical 
load. 

Recently much attention has been paid on the magnetic coupling (MC) 
process, where energy and angular momentum are transferred from a fast- 
rotating BH to its surrounding disk by virtue of the closed field fines connecting 
them (Blandford 1999; Li 2002a, hereafter L02; Wang et al. 2002b, 2003a, 
hereafter W02, W03a, respectively). 

In W02 we worked out an equivalent circuit to calculate the BZ and MC 
powers. Very recently we discussed the condition for the coexistence of the BZ 
and MC processes (CEBZMC), and found that the state of CEBZMC always 
accompanies the screw instability of the magnetic field connecting a rotating 
BH with its surrounding disk (Wang et al. 2003b, 2004, hereafter W03b and 
W04, respectively). It turns out that the screw instability will occur at some 
place far away from the inner edge of the disk, if the BH spin and the power- 
law index for the variation of the magnetic field are greater than some critical 
values. 

In this paper a new scenario for extracting rotational energy of the disk 
matter is proposed by considering the configuration of the magnetic field re- 
stricted by the screw instability. To facilitate description, this mechanism is 
referred to as the DL process, implying that energy and angular momentum 
are extracted magnetically from disk to load. By using another equivalent 
circuit we derived the expression for the DL power and torque under some 
assumptions on the unknown astrophysical load. 

This paper is organized as follows. In §2 the restriction of the screw in- 
stability to the configuration of the magnetic field is discussed, in which the 
three energy mechanisms are contained. In §3 the expressions for the powers 
and torques of the three mechanisms are derived in the two kinds of equivalent 
circuits. We compare these powers with the variation of the two parameters, 
i.e., the BH spin and the power - law index of the magnetic field on the disk. 
It turns out that the DL power is generally less than the BZ and MC powers, 
and it is comparable with the latter two when the BH spin approaches unity. 
In §4 the radiation from a quasi - steady thin disk is discussed in detail by 
applying the conservation laws of mass, energy and angular momentum to the 
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regions corresponding to the MC and DL processes. In §5 the poloidal cur- 
rent densities flowing from the BH magnetosphere into the horizon and disk 
are calculated in the regions corresponding to the three energy mechanisms. 
Finally, in §6, we summarize our main results and argue that this model can 
be used to extract clean energy from a rotating BH for powering GRBs. 

Throughout this paper the geometric units G = c = 1 are used, and the 
assumptions for BH accretion disk and the magnetic field are adopted as given 
in W03b and W04. 



2 RESTRICTION OF SCREW INSTABILITY TO CONFIGU- 
RATION OF MAGNETIC FIELD 



It is well known that the magnetic field configurations with both poloidal 
and toroidal components can be screw instable (Kadomtsev 1966; Bateman 
1978). According to the Kruskal-Shafranov criterion (Kadomtsev 1966), the 
screw instability will occur, if the toroidal magnetic field becomes so strong 
that the magnetic field line turns around itself about once. Recently some 
authors discussed the screw instability with different conditions in the BH 
magnetosphere (Gruzinov 1999; Li 2000a; Tomimatsu et el. 2001). In W04 we 
discussed the screw instability of the magnetic field in the MC process, and 
argued that the instabihty will occur if the following criterion is satisfied, 

{27rwjL)B^^/Bl<l, (1) 



where L is the poloidal length of the closed field line connecting the BH with 
the disk, and Bj^ are the poloidal and toroidal components of the magnetic 
field on the disk, respectively, and tu^ is the cylindrical radius on the disk and 
reads 



3^-6 



(2) 



In equation (2) M and a* are the BH mass and spin, respectively. The param- 
eter ^ = r/vms is the radial coordinate on the disk, which is defined in terms 
of the radius of the marginally stable orbit, and rms is related to M and a* 
by the following relation (Novikov & Thorne 1973), 
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rms = Mxi 

Xms = {3 + ^ ± [(3 - A) (3 + ^1 + 2A^f]"' , 

= 1 + (1 - alf [(1 + a,f^ + (1 - a,)^/2l 
A2^{3al + Alf' 



(3) 
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As t r o p h y s i c a I Load 




Fig. 1. Configuration of the magnetic field restricted by the screw instability 

It is found, from the criterion (1), that the screw instabihty will occur, 
provided that the BH spin a* and the power-law index n are great enough. For 
the given values of a* and n we can determine the disk region for the screw 
instability by using the criterion (1) as follows. 



^s<^<oo, 



(4) 



where = T^s/T^ms is the minimum radial coordinate for the screw instability. 
We think that the screw instability prevents the closed field lines to reach the 
region indicated by inequality (4). It seems reasonable that the configuration of 
the magnetic field in this region might consist of the open field lines connecting 
the disk with the astrophysical load, and this region is referred to as the DL 
region corresponding to the DL process. Considering the restriction of the 
screw instability, we modify the configuration of the magnetic field as shown 
in Figure 1, and the three mechanisms of extracting energy magnetically from 
the BH accretion disk are included: (1) the BZ process involved the open 
field lines connecting the BH with the astrophysical load, (2) the MC process 
involved the closed field lines connecting the BH with the surrounding disk, 
and (3) the DL process involved the open field lines connecting the disk with 
the astrophysical load. 

Now we give a brief description for the configuration in Figure 1. In W03b 
we discussed the angular boundary between the open and closed field lines on 
the BH horizon. It is shown that the boundary angle 6m exists in CEBZMC, if 
the parameters a* and n are great enough. In W04 we argued that the state of 
CEBZMC always accompanies the screw instability, and the minimum radial 
coordinate = ''^'sh'ms for the screw instability can be determined by the 
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criterion (1). Considering the restriction of the screw instabihty to the closed 
field lines and the action of the magnetic pressure on the horizon, we think 
that the boundary angle will be extended from 6m to Os- The angles Om and 
6s are related respectively to infinitive and the minimum radial coordinate 
by the mapping relation (W04): 

cos^ — cos^j:, = j G {a^;^,n) d^, (5) 



where 



G (a*; ^, n) = — . ■ (6) 

2^(1 + a^x-t + 2alx^,) (1 - 2x^1^-^ + alx^i^-^) 



Therefore the angular region of the open field lines for the BZ process is given 

by 

< ^ < ^5, (7) 



and the angular region of the closed field fines for the MC process is given by 

es<e< Ol, (8) 



where is the lower boundary angle of the closed field lines. Throughout this 
paper 6i = 0.457r is assumed in calculations. Accordingly the value range of 
the radial coordinate for the MC process is given by 

K e < ^5, (9) 



where ^ = 1 and ^5 correspond to and 9s by the mapping relation (5), 
respectively. 

Following Blandford (1976), we assume that the poloidal magnetic field 
varies with the parameter ^ on the disk as a power law ( W03a, W03b) , 

K = {rj^n {rms)\ 1< e < ^5, (10) 



where B\ is the poloidal component of the magnetic field on the horizon, 
and n is the power-law index of 5^ varying with ^. The quantities ru and 
{j'ms) are the radius of the horizon and that of the inner edge of the disk, 
respectively. 
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Fig. 2. Equivalent circuit for the BZ and MC processes (circuit I) 

According to the above discussion the value range of the radial coordinate 
for the DL process is given by equation (4), and the corresponding poloidal 
magnetic field varies with the parameter ^ as follows, 



Cs<C<oo. 



(11) 



Equation (11) is similar to equation (10), where zUj^ (r^s) and = 1 are 
replaced by and ^5, respectively. 



3 POWER OF MAGNETIC EXTRACTION FROM BH ACCRE- 
TION DISK 



It is tempting for us to discuss the power extracted from BH accretion 
disk based on the above magnetic field configuration. In W02 we derived the 
expressions for the BZ and MC powers by using an equivalent circuit as shown 
in Figure 2 (henceforth circuit I), where segments PiSi and Qi-Ri represent 
the two adjacent magnetic surfaces, and segments PiQi and RiSi represent 
the BH horizon and the loads (either the astrophysical load or the disk load) 
sandwiched by the two magnetic surfaces, respectively. 

By using circuit I we derive the expressions for the BZ and MC powers and 
those for the BZ and MC torques corresponding to the modified configuration 
of the magnetic field as follows, 

WP. = 2«;/ ,L(^_'^)^,„., , o<e<fe (12) 
Wr.^4„.(i.,)/il=i^g^, o<o<0s (13) 
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Pmc Po = 2al r'\,\. es<e<eL 

J 2 — (1 — q) sm 



(14) 



TMc/To = Aa4i + q) I }^ /^'^.^f, , es<e<e, (is) 

2 — (1 — gf) sm 



where APbz, ^Tbz, ^Pmc and AT^c are related by 



(16) 



In equations (12) — (15) the parameters k and (3 are respectively the ratios 
of the angular velocities of the open and closed field lines to the angular 
velocity of the horizon and read 



2 (1 + g) 



X 



1 -1 



+ a* 



(17) 



where Jl^ is the angular velocity of the magnetic field line, Qh and fto are 
respectively the angular velocities of the horizon and the disk and read 



2r 



n 



D 



H 



(18) 



Usually k — 0.5 is taken for the optimal BZ power (MT82), while P 
depends on the BH spin and the place where the field line penetrates on the 
disk. In equations (12) — (15) we use the parameter q = ^Jl — a^, and the 
parameters Pq and Tq are defined as 



Po = {B^Hf ^ Bl {M/Mof X 6.59 x W^erg ■ s~^, 



To = {B^H? ~ B'i {M/MqY X 3.26 x 10^^^ ■ cni^ ■ s 



(19) 



where P4 is the strength of the poloidal magnetic field on the horizon in units 
of 10^ gauss. 

Considering that the poloidal magnetic field will exert a braking torque 
on the current in the DL region, we can calculate the DL power and torque by 
using the equivalent circuit as shown in Figure 3. Henceforth this equivalent 
circuit is referred to as circuit II, where segments P2S2 and Q2R2 represent 
the two adjacent magnetic surfaces consisting of the open field lines in the 
DL region, and segments P2Q2 and R2S2 represent the disk surface in the 
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Fig. 3. Equivalent circuit for the DL process (circuit II) 

DL region and the load sandwiched by the two adjacent magnetic surfaces, 
respectively. The quantities AZa and Ae^ (the subscript "i" is omitted) are 
the resistance of the load and EMF due to the rotation of the disk, respec- 
tively. The disk load is neglected in calculations by considering the perfect 
conductivity of plasma. 

The following equations are used to derive the DL power in circuit II, 
which are similar to those given in deriving the BZ power and MC power in 
W02. 

APDL = lhL^ZA, lDL = Ae^/AZA, Ae^ = -{A^D/27i)nD, (20) 

where Idl is the current in each loop of circuit II. The minus sign in Aed arises 
from the direction of the magnetic flux between the two adjacent magnetic 
surfaces, which is given by 

A*^ = B£2n (wp/y/A)^^^^^ dr. (21) 

The concerned Kerr metric coefficients are given by (Thorne, Price & Mac- 
Donald 1986) 

w = (S/p) sin^^, 

^ ' (22) 

= cos^ 6*, 

A = - 2Mr. 



Since the astrophysical load remains unknown, we give some simplified 
assumptions as follows. 

(1) The load is axisymmetric, being located evenly in a plane P with 
some height above the disk. 
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(2) Each open field line intersects with the disk and the plane P at the 
cylindrical radii r and r', respectively. The both radii are related by 



Ar, (23) 



where A is assumed to be a constant. 

(3) The surface resistivity of the unknown load is assumed to obey the 
following relation, 

(Ji^^a^a^^ATra^, (24) 



where is a parameter, and a^j — Att — 377 ohm is the surface resistivity 
of the BH horizon. Throughout this paper ck^ = 1 is assumed in calculations. 
The load resistance AZa between the two adjacent magnetic surfaces can be 
written as 

AZA = a,^ = 2a,^, (25) 



where equations (23) and (24) are used in the last step. Incorporating equation 
(11) with equations (20) — (25), we have 

APW^o = /(a*,e,^)^^e, (26) 



where the function / (a*, ^, n) is expressed by 



Similarly, the DL power is related to the DL torque by 

APdl = ^dATdl, (28) 

and we have 

ATDL/To^g{a.,^,n)d^, (29) 



where the function g (a*, ^, n) is expressed by 
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(b) 

Fig. 4. The BZ power (solid line), the MC power (dashed line) and the DL power 
(dotted line) (a) versus the BH spin a* with n = 5.5, and (b) versus the power-law 
index n with a* = 0.9. 

Integrating equations (26) and (29) over the DL region, we have the DL power 
and torque expressed by 



PDL{a*,^,n)/Po = j f{a^,^',n)d^', 



(31) 



Tdl (a*, t n)/To j g (a*, n) . 



(32) 



is 



For the given values of a* and n the strength of Pbz/Pq, Pmc/Pq 
Pdl/Pq are compared by using equations (12), (14) and (31) as shown in 
Figures 4. 
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Prom Figure 4 we obtain the following results, 

(1) For the increasing a* with the given n, both the BZ and the MC powers 
vary non-monotonically, attaining their maxima as the BH spin approaches 
unity, while the DL power increases monotonically. 

(2) For the increasing n with the given a*, the BZ power increases mono- 
tonically, the MC power decreases monotonically, while the DL power varies 
non-monotonically, attaining its maximum value Pol ~ 0.034 with n — 6.39. 

(3) The MC power is generally greater than the BZ power, if the power- 
law index n is not very big. The DL power is generally less than the BZ and 
MC powers, and it is comparable with the latter two powers as the BH spin 
approaches unity. 



4 RADIATION FROM A MAGNETIZED ACCRETION DISK 



4-1 Radiation from the MC region 



Radiation from a relativistic steady thin disk around a Kerr BH has been 
discussed based on the three conservation laws of mass, energy and angu- 
lar momentum with the "no-torque boundary condition" by some authors 
(Novikov & Thorne 1973, Page & Thorne 1974). By combining the MC effects 
with the above conservation laws and the "no-torque boundary condition", 
the following equation of radiation from a relativistic quasi-steady thin disk 
was derived in L02, 

rptotal 771/ _i_ T? /'QQ^ 



where and Fmc are the radiation fluxes due to disk accretion and the 
MC effects, respectively, and they are expressed by 

Anr (E+ - VtoL+f J ^ ^ J \ / J ' ^ ^ 



In equation (35) Hmc is the flux of angular momentum transferred between 
the BH and the disc by the magnetic held, and and are the speciflc 
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energy and angular momentum of a particle in the disc, respectively, and they 
read (Novikov & Thorne 1973) 

E+^(l- 2x-' + a,x^') /(l - 3x-' + 2a,x-')'^', (36) 
L+ = Mx(l - 2a,x-' + alx-') /(l - 3x-' + 2a,x-')'^'- (37) 

From equations (36) and (37) we have = Ems and L"*" = Lms for = 1 with 
X = Xms- The expression for Hmc can be worked out by using the following 
relation, 

jdTMc/de) {de/dp 

diMc/or = ^ = AnrHMC: (38j 

^ A,rras 

where dTMc/dO and d9/d^ can be calculated using equations (15) and (5). In 
W03a we expressed Hmc as follows, 

Hmc /Ho = A (a„ T", K ^ < ^5, (39) 



where 

Ho = (S^)^ M = 1.48 X 10^2 X Bl (M/Mq) g ■ s'^, (40) 

A (a*, ^) = 27rx£ [2 csc^ -^^^ ('^*' > 

, (41) 



Since the magnetic field on the horizon is brought and held by its sur- 
rounding magnetized disk, there must exist some relations between the mag- 
netic field and the accretion rate. As a matter of fact these relations might be 
rather complicated, and would be very difi^crcnt in different situations. One 
of them is given by considering the balance between the pressure of the mag- 
netic field on the horizon and the ram pressure of the innermost parts of an 
accretion flow (Moderski, Sikora & Lasota 1997), i.e., 

(S^)7(87r) = Pram ~ pc' ~ M^/(47rr|), (42) 



Prom equation (42) we define Fq as 

M2 (1 + q) 



- {B^nf - . .2 - (43) 
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By using equations (34) and (35) we have the radiation fluxes, -F^^/ -Fo 
and Fmc/Fq, versus the radial parameter ^ for the given values of the power- 
law index and the BH spin as shown in Figure 5. 

Expecting Figure 5 we find that the peak of Fmc/Fq is not only closer 
to the inner edge of the disk, but also is greater than that of F^j^j Fq. Fur- 
thermore, Fmc/Fq varies more steeply with ^ than F^^^Fq does. Thus the 
radiation flux due to the MC mechanism can rcsTilt in a very steep cmissivity in 
the inner region of the disk, which is consistent with the recent XMM-Newton 
observation of the nearby bright Seyfert 1 galaxy MCG-6-30-15 (Wilms et al. 
2001; Li 2002b; W03a). 



4-2 Radiation from the DL region 



We can show that the radiation flux from the DL region also consists of 
two terms: 



where is the electromagnetic flux in the DL process, and is the 
radiation flux due to disk accretion in the DL region. The flux F^l can be 
worked out by 



Substituting equation (26) into equation (45), we have 

Fdl/Fo - (46) 



The difference between the radiation from the DL region and that from 
the MC region lies in two aspects. 

(1) The outgoing flux Fmc can be obtained after resolving the equations 
of the conservation laws with Pmc and Tmc incoming the disk, while the 
outgoing flux Fdl is given before resolving the equations from the conservation 
laws; 

(2) "No-torque boundary condition" can be used at rms for the solution 
in the MC region, while this boundary condition is not valid at for the 
solution in the DL region. 
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(c) 

Fig. 5. The radiation flux from the MC region, Fj^^/Fo (sohd line) and Fmc/Fq 
(dashed line) versus ^ for 1 < C < ^5 with n = 5.5 and different values of the BH 
spin: (a) a* = 0.5, (b) o* = 0.7 and (c) a* = 0.998. 
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Thus we find the flux by applying the conservation laws of energy 
and angular momentum to the DL region, i.e., 

^ (MnE+ - 27rrW^nn) = Anr (Fi,\E+ + HnL^o) , (47) 



dr 

d_ 
dr 



|- (m^L+ - 27rrW;) = 47rr (fI\L+ + Hol) , (48) 



where accretion rate is regarded as a constant for a quasi-steady disk, 
and is the internal viscous torque per unit circumference. Incorporating 
equations (47) and (48), we have 

- i-dnn/dr) ^^^^ 



The quantity H^l is the flux of angular momentum, which is related to the 
electromagnetic torque by 

dToL/dr = {dTDL/dO (d^/dr) = AnrHoL- (50) 



Substituting equation (32) into equation (50), we have 

Hdl/Ho = g{a.,^,n)/{47r^xL)- (51) 

Substituting equation (49) into equation (48) and integrating equation (48) 
over the DL region, we have the expression for as follows, 

FI/^^Fa + Fb + Fc, (52) 
where 

Anr {E+ - VtuL+f J ^ J \ / J ' ^ ^ 



{-dQp/dr) 



r 

J (e+ - VtDL+) rHoLdr, (54) 



-diln/dr {E+ - nnL^fvFlrl,'^' 
^ r{E+-VLDL+Y^ -dVLo/dr ^"=''^ ^ ' 

Expecting equations (52) — (55), we find the following characteristics of Fjj\: 
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(c) 

Fig. 6. The radiation flux from the DL region, Fjj\/FQ (sohd line) and Fdl/Fq 
(dashed line) versus ^ for ^5 < ^ < 4^5 with n = 5.5 and different values of the BH 
spin: (a) a* = 0.5, (b) a* = 0.7 and (c) a* = 0.998. 



16 



(1) The flux Fa has the same form as F^^ except the integral region; 

(2) The flux Fb has a similar form to Fmc except the minus sign in equa- 
tion (54), which implies the contribution of the DL torque on is negative; 

(3) The flux Fc represents the effect of the boundary condition at r = r^, 
where "no-torque boundary condition" is not valid. 

By using equations (52) and (46) wc have the radiation fluxes, F^j^/Fq 
and F]:,l/ Fq, versus the radial parameter ^ for the given values of the power- 
law index and the BH spin as shown in Figures 6. 

Inspecting Figure 6, we have the following results. 

(1) The radiation flux Fcl/Fq decreases monotonously and very steeply 
with the increasing ^. This result arises from two aspects, i.e., both the angular 
velocity decreases and the area of the ring from r to r + dr increases with the 
increasing ^. 

(2) The radiation flux F^a/ Fq generally decreases monotonously with 
the increasing ^, while it varies non-monotonously with ^ as the BH spin ap- 
proaches unity as shown in Figure 6c. This result arises from the conservation 
laws of energy and angular momentum and the behavior of Fb l/Fq near the 
boundary at r^. 



5 CURRENT DENSITIES FLOWING FROM MAGNETOSPHERE 
INTO HORIZON AND DISK 



In W03b we calculated the current density flowing from the BH mag- 
netosphere into the MC region on the horizon in the state of CEBZMC by 
using the circuit I. In this paper the three energy mechanisms of the magnetic 
extraction are described by using the circuits I and II, based on which the 
current densities flowing from the BH magnetosphere into the horizon and the 
disk can be calculated and compared. As argued in W03b the poloidal current 
flowing in each loop of circuit I can be written as 

/f.(a.^)^/o ,j4^_-f_^^ , 0<.<.„ (56) 
7^g(a.,g,n)^/o ^ es<e<eB. (57) 

where /q = ^ 1.48 x W^B^^ (M/Mq) A. The currents /f^ and I^^ are 

respectively the poloidal currents flowing on the BZ and MC regions of the 
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(b) 

Fig. 7. Current densities versus with a* = 0.9, 0l = 0.457r, (a) j§z/jo ^'^^ 
< 9 < 9s, where A, B, C are the end points of the curves corresponding to 
n=4, 5.5, 7, respectively, (b) jMc/jo iov 9s < 9 < 9l ioi n = 4, 5.5, 7 and 9 in 
dashed, sohd, dotted and dot-dashed Hues, respectively. 

horizon, which depend on the angular coordinate 6 for the given values of a* 
and n. Inspecting equations (56) and (57), we find that 1^^ ^^^^ generally not 
equal to I^q lor (3 ^ k. The conservation of the current at 9s is guaranteed 
by the current fiowing between the magnctosphere and the horizon. The con- 
tinuity of the current at the boundary between the BZ and MC regions of the 
horizon is discussed in W03b. 

Based on the conservation of current, the current densities fiowing from 
the BH magnetosphere into the above regions of the horizon are expressed by 



3bz 



• 4(l-fc)MOg cos e 
[2-smH{l-q)Y ' 



de 



O<0<0s, 



(58) 
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Fig. 8. Current densities versus 6 with n = 5.5, 9l = 0.45?!, (a) jsz/joi (b) JMc/h 
for a* =0.36 , 0.62, 0.9 and 0.998 in dashed, dotted, sohd, and dot-dashed Hues, 
respectively. 



Jmc {a*,0,n) 



27r(.z7p),,^,,^^ d0 



■JO 2-sin^ 9{l-q) [2-sin^ 0(1-?) "^^^^dg. 



, es<e<eL, 



(59) 



where = S^/(27rM) = 0.108 x B4 (M/Mq)~^ A ■ cm . For the ffiven value 
of and n, both j^^ and j^Q vary with angle 9 in the range < 9 < 9s and 
9s < 9 < 9l, respectively. 

According to circuit I the poloidal current flowing on the MC region of 
the disk is equal to that flowing on the MC region of the horizon, i.e., 

Imc n) = Igc - lo ^ ^^l; f ^_ . , 1<C<U (60) 
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Fig. 9. Current densities versus ^ with a* = 0.9, 6l = 0.457r, (a) Jmc/Jo^ (b) ioi/io 
for n =4, 5.5, 7 and 9 in dashed, sohd, dotted and dot-dashed hnes, respectively. 



Therefore the current density flowing from the BH magnetosphere into the 
disk is 



(61) 



Substituting equation (60) into equation (61), we have the current density 
flowing from the magnetosphere into the MC region of the disk as follows. 



Jmc — ^'csJe-i'i-q) 



4(1-/3) csc^gcote rfe _ d(3_ 
2csc2e-(l-Q) 



K e < 6, (62) 
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Fig. 10. Current densities versus ^ with n = 5.5, 9l = 0.457r, (a) JmcIJo^ (b) 3 dl I Jo 
for a* = 0.36 , 0.62, 0.9 and 0.998 in dashed, dotted, sohd and dot-dashed Hues, 
respectively. 



where 



1 + e 



2e 



-1^-2 



Am5 * 



(63) 



Applying the same procedure to circuit II, we can derive the current 
flowing on the DL region of the disk and the current density flowing from the 
BH magnetosphere into the disk as follows, 

(l + g)£'(a.,C,n) ^ , ^ , _ /^^.n 
Idl = -h 7. > ^5 < ^ < oo (64) 
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(65) 

^ _ ,,il+,)Cia„i) 9Dia i,n) ^ ^5 < C < OO, 



where C (a*, ^) is expressed by equation (63) and D (a*, ^) is expressed by 

^ ^ - ^S+^(«3/.^3^^+„,)Vl+S-,-ta^-2«-ix- ^^^^ 



By using equations (58) and (59) we have the current densities j^^ and 
Jmc varying with the angle 6 on the horizon for a* = 0.9 and the different 
values of n, and those for n = 5.5 and the different values of a* as shown in 
Figures 7 and 8, respectively. 

By using equations (62) and (65) we have the current densities j^Q and 
Jdl varying with the radial parameter ^ on the disk for a* = 0.9 and the 
different values of n, and those for n — 5.5 and the different values of a* as 
shown in Figures 9 and 10, respectively. 

From Figures 7 — 10, we find that the distribution features of the above 
current densities depend on both the power- law index n and the BH spin a*. 

The current densities and are always positive, while j^^^ and j^^ 
might change their signs in some cases, which correspond to the corotation 
magnetic surface (CRMS) as defined in W03b. 



6 DISCUSSION 



In this paper we discuss and compare the three powers of extracting 
energy magnetically by considering the restriction of the screw instability to 
the configuration of the magnetic field in BH magnetosphere. These powers 
are derived by using circuit I for the BZ and MC processes and circuit II for 
the DL process. It turns out that the DL power is generally less than the BZ 
and MC powers, and it is comparable with the two powers for the BH spin 
a* approaching unity. By using the conservation laws of energy and angular 
momentum we discuss the radiation fiux from a quasi-steady thin disk around 
a Kerr BH. Although the situation at the boundary between the MC and DL 
regions is not very transparent, we obtain the solution for the radiation by 
resolving the equation from the conservation laws. By using circuits I and II 
we discuss the distribution of the currents and the current densities on the 
BH horizon and on the disk for the three energy mechanisms. 
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Recently Li (2000b) presented a model for extracting clean energy from 
a Kerr BH surroimdcd by a dense plasma torus, and argued that this model 
may be relevant to GRBs, provided that the magnetic field is strong enough. 
Li argued that the baryonic contamination from the plasma in the torus is 
greatly suppressed by the magnetic confinement. In fact, the problem of the 
baryonic contamination can be also suppressed by the closed field lines in the 
MC region in our model, provided that the magnetic field is strong enough. 
Inspecting equations (12) and (19), and Figure 4, we find that the maximum 
of the BZ power could be in the following range, 

O.OSPo < PbT < 0-1^0, (67) 

which implies that P^l^ could attains ~ lO^^erg' • for M — 7Mq and 
-B4 = 10^^. Although P]^z^ seems much more than needed for powering the 
beamed GRBs, we can adjust the value of Pbz to fit the beamed GRBs by 
taking the adequate values of M and i?4 in our model. In addition, the DL 
process discussed in this paper could be another energy mechanism for extract- 
ing clean energy to GRBs. We shall apply this model to GRBs and present a 
detailed calculation in the future work. 

Acknowledgments. This work is supported by the National Natural Science 
Foundation of China under Grant Numbers 10173004, 10373006 and 10121503. 



References 

Bateman G. MHD Instabilities, 1978, (Cambridge: The MIT Press) 

Blandford R. D., 1976, MNRAS, 176, 465 

Blandford R. D., & Znajek R. L. 1977, MNRAS, 179, 433 

Blandford R. D. 1999, in Sellwood J. A., Goodman J., eds, ASP Conf. Ser. 

Vol. 160, Astrophysical Discs: An EC Summer School, Astron. Soc. Pac, 

San Francisco, p. 265 
Gruzinov A. 1999, astro-ph/9908101 
Kadomtscv B. B. 1966, Rev. Plasma Phys., 2, 153 

Lee H. K., Wijcrs, R. A. M. J., & Brown, G. E., 2000, Phys. Rep., 325 83 
Li L. -X. 2000a, ApJ, 531, Llll 
Li L. -X. 2000b, ApJ, 544, L375 
Li L. -X. 2002a, ApJ, 567, 463(L02) 
Li L. -X. 2002b, A&A, 392, 469 

MacDonald D., & Thorne K. S. 1982, MNRAS, 198, 345 (MT82) 

Moderski R., Sikora M., Lasota J. P., 1997, in Ostrowski M., Sikora M., Made- 

jski G., Belgelman M., eds, Relativistic Jets in AGNs. Uniw. Jagiellonski, 

Krakow, p. 110 



23 



Novikov I. D., & Thome, K. S., 1973, in Black Holes, 

cd. Dewitt C, (Gordon and Breach, New York) p. 345 
Page D. N., Thorne K. S., 1974, ApJ, 191, 499 
Rees M. J., 1984, ARA & A, 22, 471 

Thorne K. S., Price R. H., & Macdonald D. A. 1986, Black Holes: The Mem- 
brane Paradigm, 
Yale Univ. Press, New Haven 
Tomimatse A., Matsuoka T., & Takahashi M. 2001, Phys. Rev. D64, 123003 
Wang D.-X., Lei W.-H, & Xiao K., 2002a, ApJ, 580, 358 
Wang D.-X., Xiao K., & Lei W.-H. 2002b, MNRAS, 335, 655 (W02) 
Wang D.-X., Lei W.-H., & Ma R.-Y. 2003a, MNRAS, 342, 851 (W03a) 
Wang D.-X., Ma R.-Y., Lei W.-H., & Yao G.-Z., 2003b, ApJ, 595, 109 (W03b) 
Wang D.-X., Ma R.-Y., Lei W.-H., & Yao G.-Z., 2004, ApJ 601,1031(W04) 
Wilms J. et al. 2001, MNRAS, 328, L27 



24 



